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The  Effect  of  Accel era tion  Time  on  Plastic 
Deformation  of  Beams  under  Transverse  Impact  Loading^ 

hy 

2 

D.  St  Green 


ABSTRACT 

This  paper  considers  a special  case  of  point  impact 
loading  on  a beam  which  is  free  to  move  in  the  direction  of 
striking*  The  impact  is  such  as  to  xmpar  t to  the  mid -point  of 
the  beam  constant  acceleration  for  a certain  period,  and  then 
constant  velocity  for  all  subsequent  time.  If  this  acceleration 
is  sufficiently  great,  the  beam  will  undergo  permanent  deforma- 
tion. The  analysis  shows  how  the  magnitude  of  the  deformation 
depends  on  the  final  velocity  and  the  duration  of  the  accelera- 
tion, assuming  that  the  beam  behaves  in  a perfectly  plastic- 
rigid  manner.  This  is  the  simplest  problem  introducing  the 
finite  time  of  acceleration  that  can  be  solved,  and,  although 
it  is  highly  idealized,  its  solution  will  be  valuable  as  a 
guide  in  the  planning  of  experimental  programs. 


1.  The  results  in  this  paper  were  obtained  in  the  course  of 
research  conducted  under  Contract  N7onr-35801  betv/een  Brown 
University  and  the  Office  of  Naval  Research. 

2.  Research  Assistant,  Division  of  Applied  Mathematics,  Brown 
University. 
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List  of  symbols 

m = mass/unit  length  of  beam 

2.K,  = length  of  beam 

MQ  = fully  plastic  moment  of  beam 

x = distance  along  beam  from  mid-point 

t = time 

t = duration  of  acceleration  pulse 

ts  = time  at  which  outer  hinges  disappear 
T = total  time  required  to  compiece  deformation 
vQ  = velocity  of  mid-point  of  beam 

v = velocity  of  point  distant  x from  mid-point  of  beam 
V = final  velocity  of  beam 

R = striking  force  applied  at  mid-point  of  beam 
M = bending  moment  at  x 

\ = distance  of  outer  hinge  from  mid-point  of  beam 

K = Wl 

= initial  value  of 
K = (1  - 20k 

9 = angle  of  rotation  of  inner  segment 

0 = final  magnitude  of  9;  when  t = T,  9-0 

q>  = angle  of  rotation  of  outer  segment 

$ = final  magnitude  of  tp. 
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INTRODUCTION 

The  assumption  of  plastic-rigid  behavior  has  been 
adopted  several  times  in  previous  analyses  of  the  plastic  deform- 
ation of  beams  subjected  to  various  types  of  impact  loading 
[1  - 4]^.  It  is  now  desirable  to  verify  its  practical  applica- 
bility by  experiment.  This  can  only  be  done  if  the  theoretical 
solutions  relate  to  experimental  conditions.  In  impact  tests  it 
is  easier,  for  example,  to  measure  velocities  than  forces.  With 
this  in  mind,  Symonds  and  Leth  [l]  solved  the  proW  em  in  which 
the  mid-point  of  a free  beam  is  suddenly  caused  to  move  with  a 
given  velocity;  they  neglected,  however,  the  effect  of  finite 
time  of  impact.  Our  present  analysis,  of  which  theirs  is  the 
special  case  of  zero  acceleration  time,  introduces  this  factor 
and  therefore  approaches  more  closely  to  real  conditions.  It  is 
our  intention  here  to  provide  information  on  the  effects  of 
acceleration  and  pulse  time  which  will  be  useful  in  the  arrange- 
ment and  interpretation  of  tests. 

The  problem  solved  in  this  paper  will  now  be  specified. 
A uniform  beam  AB  of  length  2 £ is  at  rest  and  without  lateral 
restraint,  when  its  mid-point  0 is  struck  laterally  in  such  a 
way  as  to  give  to  0 constant  acceleration  for  time  m,  and  con- 
stant velocity  V thereafter  (Fig.  1).  It  is  required  to  find  the 
resulting  angle  of  deformation  at  the  mid-point  of  the  beam 
(0  in  Fig,  2),  assuming  that  the  relation  between  bending  moment 
and  curvature  is  that  shown  in  Fig.  3.  Perfectly  plastic-rigid 

3^  Numbers  in  square  brackets  refer  to  the  Bibliography  at  the 
end  of  the  paper. 
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me chani  cal  behavior  means  that  there  are  no  elastic  effects  and 
the  beam  bends  only  where  it  is  subjected  to  the  fully  plastic 
moment  MQ,  at  which  sections  it  can  assume  unlimited  curvature 
while  M0  is  maintained*  Thus  elastic  vibrations  are  precluded 
and  the  solution  is  conspicuously  simplified.  An  analysis  based 
on  this  assumption  is  likely  to  produce  valid  results  for  ductile 
materials  only  when  the  total  energy  absorbed  in  plastic  deforma- 
tion exceeds  by  far  the  total  elastic  energy  the  beam  could  store 
in  bending  [1,2], 

There  are  three  possible  types  of  configuration  that 
cou3.d  arise  from  the  motion  imparted  to  the  beam.  Which  of 
these  obtains  depends  as  follows  on  the  magnitude  of  the  non- 

p 

dimensional  quantity  mt  V/Mdt  : 

p 

(a)  m-t  V/M0t  < 2,  where  m represents  the  mass  per  unit  length  of 
the  beam.  The  limit  moment  M0  is  never  achieved  and  hence 
there  is  no  deformation. 

(b)  2 < mJt?V/K0%  _<  41.7  A plastic  hinge  appears  at  the  mid- 
point and  all  the  deformation  is  concentrated  there. 

(c)  m-L  v/M0t  > 41.7  Three  plastic  hinges  occur,  one  at  the 
mid-point  and  two  symmetrically  placed  outer  ones.  The 
initial  positions  of  the  outer  hinges  depend  upon  V/t,  but 
these  hinges  move  outwards  until  in  the  course  of  time  they 
r-ach  points  a quarter  of  the  length  of  the  beam  from  the 
ends,  and  then  the  bending  moment  is  no  longer  sufficient  to 
sustain  mhem  and  they  disappear  (Fig,  'h ).  These  outer  hinges 
leave  behind  them  smoothly-deformed  regions* 


Conditions  are  always  symmetrical  about  the  mid-point  of  the 
beam.  After  stating  briefly  the  assumptions  made  to  facilitate 
our  analysis,  we  shall  discuss  cases  (b)  and  (c)  in  detail.  Case 
(a)  is  trivial. 

Apart  from  the  neglect  of  elastic  effects  the  principal 
assumptions  made  are  that  the  limit  moment  is  unaffected  by 
strain-rate  [5],  strain-hardening  and  shear  stresses  [6,  ?]. 
Experiments  are  needed  to  tell  us  how  important  these  phenomena 
are,  and  within  what  range  our  results  are  useful.  It  is  perhaps 
fortunate  that  although  strain-hardening  and  high  strain-rate 
tend  to  increase  the  limit  moment,  shear  stresses  influence  it  in 
the  contrary  sense,  A further  assumption  made  throughout  the 
work  is  that  the  changes  of  shape  of  the  beam  are  so  small  as 
to  have  a negligible  effect  on  the  equations  of  motion. 


ANALYSIS 

1.  Case  of  One  Plastic  Hinge  (2  < m-£2V/M0T  _<  ^1,7) 

2 

If  V/t  is  greater  than  2M0/m-t  the  limit  moment  is 
always  reached  at  the  mid-point  of  the  beam  and  a plastic  hinge 
occurs  there.  This  can  easily  be  shown  by  considering  the  bending 
moment  distribution  caused  by  inertia  forces  in  a rigid  rod. 

The  equations  for  the  rotation  of  either  half  of  the 
beam  (Fig,  5)  during  the  pulse  and  for  the  subsequent  time  are 
respectively 


and 


•• 


© = 


3 v. 

7 * " raJ 


(1) 


© = - 


3M 
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m-t 


t 


(2) 
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The  plastic  hinge  disappears  when  9=0.  Application  of  this 
condition  to  (1)  and  (2)  shows  that  the  final  angle  of  deforma 
tion  developed  by  the  beam  is 


© 


M 'v 
o 

miV 


(3) 


p 

The  value  of  mt  V/Mqt  at  which  throe  plastic  hinges 
first  appear  can  be  elicited  from  the  bending  moment  distribution. 
Let  M be  the  bending  moment  at  a point  distant  x from  0.  The 
equation  for  the  translatory  motion  of  the  beam  during  the 
acceleration  pulse  is 


whence, using  (1), 


R 

7 


- © 


R 

2 


(4) 

(5) 


where  R is  the  total  force  applied  to  the  mid-point.  The  moment 
at  0 is  the  limit  moment  M0,  therefore,  reckoning  sagging  bend- 
ing moment  positive, 

x 

P 

M = Rx/2  - Mc  - m i (x  - x*)(V/t  ••  x*  ©)dx*,  (6) 

Jo 

• • 

By  substituting  for  R and  9 and  then  finding  the  value  of  x at 
which  M is  a maximum,  wo  find  that  the  outer  hinge  first  appears 
when  mt2V/M0T  = 41.7  and  then  x/Z  = 0.404, 

In  Table  I are  shown  some  values  of  the  final  angle  0 
for  corresponding  values  of  m-t2V/M0a-  lying  between  2 and  41.7* 
These  have  been  plotted  in  Fig,  9* 
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TABLE  I 


n 

mt2V  ; 

M % 
o 

! 

i 

2 ; 3 

! 

— - -j  - 

i 

H 

o 

20 

30  j 41,7 

l 

M © i 

0 I 

m-tV2 

L 

1 

o ; 0.125 
! 

CO 

00 

rH 

. 

o 

; 

i 

0.225  : 0.300 

! 

* 

0.338 

u.. 

0.351  1 0.357 

! 

j 

2 

2,  Case  of  Three  Plastic  Hinges  (m-l  V./M0t  > 4l6?) 

If  M/%  is  greater  than  41.7  M three  plastic  hinges 

form  simultaneously  at  the  common cement  of  the  motion.  The 
initial  positions  of  the  symmetrically  situated  outer  ones  depend 

on  the  magnitude  of  M/% . We  shall  now  seek  equations  of  motion 

• • 

for  the  beam  by  first  taking  any  function  ~ vQ(t)  for  the 
acceleration  of  the  mid-point.  These  general  equations  can  then 
be  suitably  adapted  to  represent  each  phase  of  the  motion  in 
turn  (e.g.  by  putting  vQ  = M/%  = constant  for  0 < t £ t). 

Conditions  being  symmetrical  about  the  point  0,  we 
need  consider  only  half  of  the  beam,  OA  (Fig.  6).  H represents 
the  outer  hinge  in  this  half.  The  general  equations  of  motion 
for  OH  and  HA  respectively  are 


2M0  = 


M - 
J1o 


rnX,c'  Vo-  ^ 9 

~ v0  -y J 

£li&  - X,)3  " 

12  Y 


(7) 

(8) 


where  \ = \(t)  = distance  of  the  outer  hinge  from  the  mid-point  0 
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© = angle  of  rotation  of  the  inner  segment  OH; 

cp  = angle  of  rotation  of  the  cuter  segment  HA, 

The  element  of  mass  mdX,  which  is  transferred  from  the  outer  to 
the  inner  segment  as  the  yield  hinge  traverses  a distance  dX 
not  influence  the  equations  because  its  momentum  is  not 
affected  by  the  transfer.  Equation  (8)  follows  from  there  being 
zero  shear  force  at  H,  and  so  from  the  acceleration  of  the  centre 

of  mass  of  HA  being  zero.  Further  equations  represent  the  veloc- 

ity and  acceleration  distribution  along  the  rod:  for  x > X 

> ? 

i • * 

v = vQ  - i © dx  - ; cpdx,  (9) 

JO  d X 

and  therefore,  by  the  rule  for  differentiation  under  integral 


signs 


v — Vq  — j © dx  - X© 


max  + Xcp , 


(10) 


But  the  acceleration  of  the  mid-point  of  HA,  x - -(-L  + X)  is 
zero,  whence 


vQ  - X©  - - \)<p  + X(cp  - 6)  = 0, 


(ID 


Elimination  of  © and  cp  from  (11)  gives 


m(l  - X)2 


= - X( cp  - 9) 


(12) 


Having  derived  the  general  equations  (7)  to  (12),  let 
us  investigate  the  motion  and  deformation  of  the  beam  in  the 
following  three  phases: 
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I during  the  acceleration  pulse  (0  < t < t)j 
II  the  time  between  the  end  of  the  pulse  and  the  elimination 
of  the  outer  plastic  hinges  H(t  < t _<  ts)  > 

III  the  time  after  this  and  until  no  more  bendjng  occurs  at  the 
mid-point  0(ts  < t .<  T), 

It  will  be  evident  from  our  subsequent  analysis  that  the  outer 
hinges  cease  to  be  sustained  before  rotation  ceases  at  the  central 
nri  p s 

Phase  IT  0 < t £ T 

In  this  phase  v = V/t  = constant.  It  can  be  shown 
that  in  this  circumstance  (12)  demands  that  X = 0 and  so  X - X^ 

- constant  for  t < t.  Thus  from  (12),  and  putting  \/£  = E • 


12 M t ^(1  - E )2 

o _ o o 

m^2V  (1  - 2Kq) 


(13) 


We  have,  therefore,  that  throughout  the  pulse  period  the  hinges 
do  not  move  along  the  rod;  the  greater  the  value  of  V/t  the 
closer  are  they  to  the  mid-point;  their  initial  positions  can 
never  be  farther  from  the  mid-point  than  O.koWL  because  V/t 

r\ 

must  exceed  ‘fl,7M0/m'Cr . Figure  7 shows  V/t  as  a function  of  £0; 
the  curve  is  asymptotic  to  the  line  K0  - 0.  Henceforth  in  our 
work  we  shall  express  deformation  as  a function  of  F instead  of 
V/t  because  of  the  consequent  simplification.  Figure  7 is  a 
means  of  rapid  transformation  into  terms  of  V/t  if  such  is 
required. 

Returning  to  (7)  and  (8),  for  OH  we  have 
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XoT  m\3 

2(1  - 2g0)  - i\ 

2r0(!  - 2t,  Q) 


(using  (13)) 


(iVa) 


(llrt>) 


1 9 x 

2 


(15) 


and  similarly  for  HA 


12Mqt 


ia->.„)3  a - 250)d  - o 


(16) 


9-  = 1 


2 r 


(17) 


Substitution  of  values  for  in  these  equations  gives  the  re- 
sults representing  conditions  at  t = % in  the  first  two  lines  of 
Table  II,  page  16.  It  is  our  object  in  this  paper  to  determine 
the  resultant  deformation  ©,  so  only  the  figures  relating  to  it 
have  been  tabulated,  and  those  for  9 have  been  omitted. 

Phase  II,  t < t < 

The  deflections  and  velocities  at  t = t are  the  initial 
conditions  for  the  differential  equations  of  motion  for  this 

phase.  The  acceleration  condition  now  is  that  v =0,  and  so 

o 7 

(7),  (8)  and  (11)  become  (18),  (19)  and  (20): 


9 = - 


— & -L 

m-L3  £3 


m£3  (1  - O' 


(18) 


(19) 


- X©  - SA  r-.M  9 + X(  9 - 9)  = 0 


(20) 
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Elimination  of  9 and  9 from  those  equations  gives 


i - 0 = j C -L—  - 4_  ] 

^ (l  - o v 


(21) 


Differentiation  of  (21)  with  respect  to  t,  and  then  a further 

• « • • 

suVstitution  for  9 and  9 yields  the  following  equation  relating 


( to  t: 


i _ (l  - o£ 
r-  2 ~ Tl 2t  - : 


(22) 


the  solution  to  wnich  is 


KK/la  - 20  exp( £/2 ) ] =--  A 


(2j) 


where  A is  a constant  of  integration  depending  on  the  angular 
velocities  of  the  segments  at  t = t, 


72M  C 

A = 

m£2V 


(24) 


C=CU0)  = (l-2$0)  /6(l-^0)[3(l-250)(l-50)-2^-6(l-C0)-"]exp(50/2) 

(2?) 

The  relation  C = C(qQ)  is  shown  in  Fig,  8.  This  phase  of  the 
problem  is  virtually  solved  now  that  (23)  represents  the  movement 
of  the  outer  hinges.  It  merely  remains  to  find  the  increase  in 

o 

© during  the  period  % < t < ts>  and  the  angular  velocities  9 and 
© at  t by  integrating  (18)  and  (19), 

From  (18)  and  (23), 


Sdt  = - 


6Mo  dt  _ 2 V I r -1  -.4.-2  2 t,.  „4xn,v 

— 7 T " 7 7?  \ (1_K  ) K exp[-  f(l-L  ) JdK 

mV  V 3 ^ ' 4 

'JK0  (26) 


l 


All -112 


-»  O 


where  K = (1  - 2£ ) and  KQ  corresponds  to  £ . Then 


Qdt  = | X “ (n  + l)K(i+n+2)exp[-  kl  - K4)]dK  (27) 

'tG  n-~-0  4 


The  fact  that  when  t = t , £ = ^ and  K = 0 follows  from  (21)  by 

s 2 

• • 

putting  cp  = Q,  which  is  the  condition  for  the  disappearance  of 
the  outer  hinges.  Further,  for  the  increase  in  9 during  this 
phase  we  have 


Qdt  = ; (0  - 0 )dt  + Q dt 


(28) 


which  has  to  be  evaluated  by  the  same  means  as  was  (26).  The 
computed  results  are  in  Table  II.  The  subscript  s refers  to 
conditions  at  t = t . 

o 

Phase  III,  ts  < t < T 

Only  the  hinge  at  the  mid-point  now  remains,  and  the 
two  halves  of  the  beam  rotate  as  rigid  bent  rods.  The  equation 
of  motion  for  0A  is  therefore 


0 = - 3M0/nvtJ 


(29) 


So  the  angular  acceleration  is  constant.  The  rotation  persists 
until  the  angular  velocity  is  reduced  to  zero.  Then  t = T and 
9 reaches  its  final  value  ©.The  results  for  this  phase  are  also 
in  Table  II. 
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Resu?.tant  Deformation  for  Whole  Period  0 < t < T 

Collecting  together  our  data,  we  have  ©,  the  final 
angle  at  the  mid-point  (when  three  hinges  have  been  present), 
given  in  Table  II.  All  these  results  are  subject  to  the  angles 
being  small,  otherwise  the  equations  of  motion  would  be  inva] i- 

dated  by  changes  of  shape.  The  working  has  been  checked  by 

« • 

graphical  integration  to  ensure  that  = Qs,  and  also  by  drawing 

up  an  energy  balance  for  the  case  of  = 0.3*  Only  0 has  been 

described  here.  Calculations  show  that  0 - $ is  much  less  than 
0.  Figure  2 indicates  qualitatively  the  sort  of  deformation  that 
occurs.  It  should  be  observed  that  most  of  the  angle  © is  built 
up  during  the  period  t < t < T. 

o * 

3 • Results 

All  the  values  of  G In  Tables  I and  II  (i:e.  for  both 
types  of  deformation  configuration)  have  been  plotted  in  Fig.  9> 
-which  contains  the  most  important  numerical  results  of  the 
analysis.  In  the  graph  the  line  M^/m-Cv2  = O.h-25  has  been  ob- 
tained from  the  work  of  Symonds  and  Leth  [1]  in  which  they 
considered  the  limiting  case  of  infinite  acceleration.  This 
result  can  also  be  derived  by  the  method  used  in  the  present 
paper:  (13)  shows  that  for  V/t — >00,  £0  — >0,  and  then  suitable 
modification  of  (23)  gives  the  required  answer.  It  is  a satisfy- 
ing check  on  the  computation  for  Fig.  9,  that  the  graph  of 

is  asymptotic  to  this  limiting  line.  It  is  worthy  of 
note  that,  provided  three  hinges  occur,  our  curvo  lies  at  all 
places  within  15%  of  the  limiting  case  and  within  this  degree  of 
approximation  it  is  independent  of  V/t. 
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4.  Illustrative  Examples 

(a)  For  a •£"  square  steel  bar  6"  long,  with  mid-section 
accelerated  to  50  ft/sec  in  1/20  millisec 

m^y/M0T  = 42 

so  v;e  have  the  three  hinge  type  of  deformation.  Figure 
9 gives 

0 = 0.358  rmtV2/ MQ. 

Therefore 

6 r 0,15  rdns. 

(Notes  the  limiting  case  of  Symcnds  and  Leth  would 
give  0 = 0.18  rdns). 

(b)  An  8’1  x 4g"  @ 23  lbs,  American  Standard  I beam 
32  ft.  long,  with  mid-section  accelerated  to  50  ft/sec 
in  1 millisec,  also  deforms  with  three  plastic  hinges 
because 

m-t2V/M0r  = 210 
Figure  9 gives 

0 = 0,395  m-tv2 /Mq 
= 0,22  rdns, 

(Note:  the  limiting  case  of  Symonds  and  Leth  would 
give  0 = 0.24  rdns.) 


All-112 


15 


CONCLUSIONS 

Figure  9 contains  the  major  results  of  this  work.  The 
type  of  deformation  that  occurs  in  a beam  subjected  to  the  load- 
ing specified  in  this  problem  is  governed  by  the  non-dimensional 

2 2 

quantity  m-V v/M^,  For  m^V/M^  _<  2 no  plastic  deformation  is 

p 

caused.  If  2 < rmi  V/Mct  < 4l.7  one  hinge  occurs  at  the  mid- 
point and  the  damage  can  be  ascertained  for  any  particular  case 

2 

by  using  Fig.  9.  When,  however,  m-l  V/Mdt  > 4-1.7  the  quantity 
Mo0/m-vV^  is  substantially  constant  and  an  approximation  to  within 
10$  is  that 

0 = 0.4  m£v2/Mo  ml2V/MDT  >41.7  (30) 

Thus  © may  be  considered  to  be  independent  of  % in  this  region 
and  to  depend  only  on  the  square  of  the  velocity  V. 

Experiments  are  now  required  for  comparison  with  the 
results  of  this  paper.  It  will  then  be  known  how  influential 
are  the  various  effects  neglected  in  the  analysis,  chief  among 
which  are  elastic  vibrations,  shear  forces,  strain-hardening, 


and  strain-rate 
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*Low  accuracy  because  of  slow  convergence  of  series  expansions  in  (27)  and  (28) 
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